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ABSTRACT 


An aeromagnetic survey in the Katmai National Monument was carried 
out by the Geophysical Institute, University of Alaska, in the summers 
of 1970-1971. Geologic and magnetic studies of the Katmai area have fur- 
ther demonstrated the close relationship between the Katmai Caldera, No- 
varupta plug, and the pyroclastic flows in the Valley of Ten Thousand 
Smokes . 

The magnetic fields observed appear to be associated with the 
thickness of the pyroclastic flow and the different rock units within it 
for lower flight levels, and also the contrast between the valley fill 
and the rock units at the Valley margins . 

Consistent magnetic anomalies are associated with the larger fuma- 
role lines, which were presumably sites of large scale activity, while 
the smaller fumaroles are not usually seen in the aeromagnetic map, A 
possible correlation between low positive anomalies and nuee ardente 
deposits was revealed by the aeromagnetic survey, but was not strong, 

Novarupta plug is represented by a large positive anomaly which is 
thought to represent both the lava dome and concealed intrusive material 
beneath it. In contrast, the rim around the plug has negative anomalies 
similar to those associated with Naknek Sediments. 

Around Novarupta, radiating lines indicate local fault zones, which 
in turn have controlled the fumarolic activity. The flank terraces also 
seem to be associated with local faulting. 

Volcanic trends in the Alaska Teninsula are thought to be controlled 
by the tectonics , of the Aleutian island arc system. The magnetic signa- 
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tures of this tectonic trend are not strong, but some features can be 
recognized from the distribution of positive anomalies at the 1220 m 
flight level. 

A ground survey was also carried out in several parts of the Valley 
with a view to detailed delineation of the magnetic signatures of the 
pyroclastic flow, as an aid to interpreting the aeromagnetic date. 
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CHAPTER I 


INTRODUCTION 

Katmai National Monument, located on the Alaska Peninsula, is part 
of the active Circum Pacific Volcanic and Earthquake zone. It is also 
part of the island arc system that extends for more than 3000 km from 
Kamchatka Peninsula to the Alaska Range. Numerous volcanic and non-volca- 
nic peaks are grouped around Mt. Katmai. The higher peaks, up to about 
2140 m are volcanic summits superimposed on the Aleutian range, which it- 
self reaches 1200 m ^ 1530 m and is covered by recent glaciers. The vol- 
canism of the Aleutian chain is typical of the Pacific Rim Andesite 
line". The volcanics range in age from Tertiary to recent, and consist 
mainly of basalts, andesites, and rhyolites. 

There are mostly andesites and dacites in the Katmai area, which 
have formed typical s tratovolcanoes and lava domes, and, in the case of 
Mt. Katmai itself, a caldera. The 1912 Mt. Katmai caldera- forming erup- 
tion was one of the largest events in historical times, and although many 
scientists have investigated the. origin and nature of the activity, there 

are still many unsolved problems (Fenner 1922, Griggs 1912, Williams 1954, 

v 

Keller and Reiser 1959, and Wilcox 1959). 

Recently, several petrological, tectonic, and geophysical studies 
have been made. Forbes et al . (1968), Ray (1967), Kienle (1967), and 
Burk (1965) . 

Work is currently being undertaken in the Mt. Katmai area by members 
of the Geophysical Institute, University of Alaska. The work described in 
this dissertation was part of this project, undertaken in the summer of 
1970. The main objects of the overall study were to study: 1) Current 
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volcanic activity ; 2) The general geology of the Mt. Katmai area, in- 
cluding reconstruction of the eruptive events of 1912, their origin, and 
the character of the pyroclastic flow in the Valley of Ten Thousand 
Smokes; 3) The thickness of the pyroclastic flow in the Valley of Ten 
Thousand Smokes using seismic, gravity and magnetic techniques; 4) The 
geomagnetic field anomalies associated with the volcanic activity; 5) The 
geomorphology and glacial features of the area* 

This thesis is concerned principally with the aeromagnetic data. 

From the magnetic measurements, iso-intensity maps and anomally maps have 
been compiled. 

The base camp was located on the flank of Baked Mountain, which is 
situated about 7 miles northwest of Mt. Katmai, and about 30 miles from 
Brooks Camp on Naknek Lake (Fig. 2). It was noted during the course of 
the investigation that the Alaska Peninsula is a storm lashed area with 
very unpredictable weather; the wind at times exceeding 50 m/sec. The 
weather is characteristically foggy, cloudy, and windy, the violence of 
wind being sufficient to deposit eolian sediments in the lower valley. 
Although the Alaska Peninsula is narrow, the Katmai National Monument 
area is a transition zone between the maritime climate of the Aleutian 
region and semi-continental climate of the terrestrial areas of Alaska. 

Geo logy of Mt . Katmai Are a 

The Alaska Peninsula consists in part of thick sequences of volcanic 
and marine geo-synclinal sediments of Permain to Early Jurrassic age, 
which have been intruded by Early Jurassic granitic plutons . Minor defor- 
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Fig. 2 


Geological map of the Katmai area 

Jn: Upper Jurassic, Naknek formation 

QTv: Quaternary (?) volcanic rocks, chief! 

Qt: Quaternary tuff, deposited during ear 



andesite. 

■ phase of 1912 Mount Katmai eruption. 



illation and uplift occurred throughout Cretaceous, and Tertiary times, 
and has continued up to the present. Most of the prominent structural 
features of the Alaska Peninsula were formed by post Miocene deformation 
(Burk 1965): 

It is said that the island arc volcanism began in the late Tertiary 
(Kagami et al 1965, Von Huene and Shor 1969) and that the main crustal 
movement occurred in late Pliocene times. 

On January 6, 1912, the peak of Mt. Katraai disappeared accompanied 
by violent activity in the area, including the eruption of material from 
vents and/or fissures during a comparatively short time (Griggs 1920). The 
erupted material consisted of pyroclastics in fragmental condition as ash, 
lapilli, and bombs, and also Nuee ardentes. Following these events, the 
National Geographic Society and U.S.G.S. sponsored a comprehensive survey 
of the area (Fenner 1922, Griggs 1923). 

Before the 1912 eruption, the only geological description near Katmai 
was that reported by J. E. Spurr in 1900. He traversed a deep valley 
(assumed to be the ancestral VaiJ.ey of Ten Thousand Smokes) floored with 
stream and glacial debris, and says that its walls consisted of marine 
sedimentary rocks containing Jurrassic fossils. These rocks in the Alaska 
Peninsula are thought to have been lifted and intruded by granitic plu- 
tonic rocks before the Cenozoic volcanoes began to erupt. 

Generally the volcanic history of the Alaska Peninsula has been 
characterized by a uniformly developing series, such as strato volcano 
erosion and/or caldera formation, followed by formation of a new central 
cone. In the Katmai case, the building of a new central cone has not yet 
occurred, perhaps because Novarupta and Trident are being formed instead. 
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The most noteable activity today is associated with Mts . Griggs, 
Martin, Mageik, Trident, and Novarupta, and there is no sign of current 
activity in the Katmai caldera. 

Mt. Griggs, Mt. Mageik, and Mt. Martin are s tratovolcanoes and in- 
clude andesitic and dacitic lavas of Quaternary age. The lavas are of 
different colors and have different mineral compositions , but are very 
similar in chemical composition. This probably indicates that they were 
derived from the same sources. In considering Novarupta, Williams (1954) 
said that ’’the intimately mingled rhyolitic and andesitic ejecta of which 
the avalanche deposit consists were formed, not by assimilation of old 
andesitic rocks in fresh rhyolitic magma, but by simultaneous discharge 
of rhyolitic and andesitic magmas from the same or closely adjacent fis- 
sures . The lithological fragments of andesite and sedimentary rocks in- 
cluded among the ejecta, and the banded lava in the dome of Novarupta 
only showed a trivial amount of assimilation”. On this basis, it would 
seem that there are two kinds of magma which were erupted together. Berg 
and Kubota (1967) show evidence that supports two magma chambers at dif- 
ferent depths in the vicinity of Baked Mountain base camp. 

In contrast to this, Fenner (1925) suggested that before and during 
the Katmai eruptions, a "lake” of rhyolitic lava dissolved approximately 
two cubic miles of the andesitic and sedimentary rocks of the flanks of 
Katmai before being erupted as a hybrid pumice. This concept was based on 
the idea of volcanic cycles (Bowen 1928) that change from basic to acidic. 

At or about the same time as the glowing avalanches were discharged, 
presumably from Novarupta and adjacent fissures at the head of the Valley 
of Ten Thousand Smokes, the entire top of Mt. Katmai volcano collapsed. 
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leaving in place of its pointed peak an oval crater almost 5,0 km long 
and 3.2 kra wide. It seems probable that a connection between Katmai and 
Novarupta existed at some depth, perhaps in the form of lava conduits. 
During the final phase of the eruption, rhyolitic lava issued sluggishly 
to build a dome within the Novarupta Crater. According to Williams, the 
main eruption of Novarupta certainly occurred after the collapse of Mt, 
Katmai. 

Williams classifies volcanic calderas as of Krakatoan type, when 
the caldera is associated with voluminous explosive eruptions of siliceous 
magmas, and of Katmai type when the collapse results from drainage of the 
central conduit of a volcano and perhaps also of some of the underlying 
reservoir by discharge of magma through adjacent conduits. He claims that 
"the top of Mt. Katmai collapsed in 1912, when its central conduit was 
drained by eruption of ash and pumice flows from vents in the adjacent 
Valley of Ten Thousand Smokes". 

Naknek Formation 

The name Naknek Formation is applied to a sequence of rocks of late 
Jurrassic age, locally having a basal conglomerate member. This basal 
member is identified by Mather (1963) as the Chisik conglomerate, based 
on the lithologic similarities of the conglomerate in Mt. Katmai area 
to the one described by Martin and Katz (1912) in the Iliamna Lake region. 
The upper part of the Naknek formation in the Katmai area consists of 

\ 

marine sandstone, silts, shales, and occasional lenticular fossil beds. 
Rocks of the Naknek formation are wide spread on the Alaska Peninsula as 
well as in other areas . The upper parts of Naknek Formation are ICimmerid- 
gian to basal Portlandian in age (See Fig. 3). 
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Tertiary Deposits 


Tertiary coal bearing rocks are present on the Kenai Peninsula, and 
east of Mt. ■ Katina i areas. They locally interfinger with volcanic rocks 
which may be of Eocene age (Keller and Kieser 1959) . The formation con- 
sists of abundant nonmarine conglomerates, sandstones, silts and coal 
bearing mudstone. These appear to be lagoonal or swampy sediments and 
were probably deposited in Eocene times. These coal bearing formations 
may be comparable with the Kenai Formation which has been disturbed by 
Miocene to recent movement. However, directly comparable rocks could not 
be found in the survey area. 

Volcanism began at least as early as Eocene time and continued through 
the Tertiary and Quaternary period to the present. Tectonic activity along 
parts of the circum Pacific belt has continued at least since Mesozoic 
times, but the main movement may have occurred between Pliocene and 
Pleistocene times. There are, however, no signs of a marked change in 
relative motions similar to those that produced the predominately Miocene 
volcanism and 11 green tuff movements 1 ' of Japan (Minato, Gorai and Funabashi 
1965). 

Surficial deposits in the Katmai area are predominately alluvial 
gravel and glacial moraine, these deposits produced the periglacial and 
proglacial areas in late glacial ages. Iluk Arm Lake and Brooks Lake show 
excellent examples of moraine deposits. 

Very recent deposits are dominated by the emplacement of the thick 
tuff deposits of the eruption of Mt, Novarupta. This has been discussed 
in detail by Griggs (1922), Fenner (1923), and Forbes (1970). 
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Tectonic and Volcanic Setting 


The eastern Aleutian Trench lies parallel to the Alaska Peninsula 
and seems to disappear beneath the Alaskan Mainland. 

The Eastern Aleutians have an arcuate shape of ^ 1380 km radius, 
centered off Nome; the trench is an asymmetric V-shaped and generally 
broad depression. It has half graben fault structure on the continental 
slope side, and is filled with flat -lying undeformed sediment (VonHuene 
and Shor 1969) . The Alaska Peninsula has three main echelon fault blocks 
which are bounded on one side by a steep reverse fault; these have formed 
the complex faulted anticlines of Pliocene age (Burk 1965) (Tig. 4) . 

Burk (1965) in his description of the geology of Alaska Peninsula 
and the adjacent continental shelf, noted that the present eastern Aleu- 
tian Trench could not be older than early Tertiary and possibly much 
younger. 

Hamilton (1966) has estimated the age of the trench as Cretaceous 
to Miocene, favoring Eocene. However, the eastern Aleutian trench is not 
very deep, so that it appears more as a trough than a trench as opposed 
to the western Aleutian island arc which has a well-developed deep trench. 

The Alaska Peninsula volcanism and associated trench have only 
shallow earthquakes indicating a low angle Benioff Zone and also thin 
undeformed sediments in the trench. This may indicate a young age. 

The distance between the trench and volcanic front is about 150 ^ 

200 km, the depth of the earthquake zone beneath the volcano line is 80 
km so that the angle of the earthquake zone can he calculated to be about 
20 degrees. 

Kanomori (1970) maintains that because of the relatively shallow 
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Late Mesozoic tectonic elements of southern Alaska, modified from 
Gate and Cryc ( 1963 ), Miller and others ( 1959 ), and Payne ( 1955 ). 
Areas of deposition are hachured, and areas apparently eroded are 
blank. Major faults are shown by heavy lines (after Burk 1965 ). 



depth, the temperature and pressure are not enough for large scale melt- 
ing of magma, hence large stresses can be built up resulting in larger 
earthquakes in the area. 

The origin of andesites from calc-alkali rock series and their bear- 
ing on island arc structures was discussed by Kuno (1963). He discusses 
the possibility of producing the successive magmas by fractionation of a 
basaltic parent magma. For example, northeast Japan may be divided into 
zones of tholeiite, high alumina basalt and alkali olivine basalt arranged 
in order from the Pacific side of the islands toward the continental side. 
In general, it is found that the pigeonitic rock series, the high alumina 
basalt series and alkali rock series occur in their respective zones. 

Using Kuno’s classification of magma, the Aleutian volcanic rocks 
were generated from a tholeiitic magma. However, the rocks belong to high 
alumina basalt (Ray 1967) , (Sugimura 1967) , so parent magma would 
have to have high alumina and alkali basalt affinities (Ray 1967) . Such 
melts could be derived from the partial fusion of a basaltic layer at 
the base of the crust. Kuno pointed out that the type of volcanic rock 
correlates with the depth of the seismic zone. If magma is generated at 
the seismic zone, tholeiitic basaltic magma would be generated at 130 ^ 

160 km, high alumina basalt at 160 'v 250 km, and alkali basalt magma at 
more than 250 km. 

On the other hand, Sugimura (1963) proposed a classification of is- 
land arc rock series using an index which is computed according to the 
formula: 

0 = Si0 2 - 47 (Na 2 0 + ^0) / A1 0 3 

where Si0 2 is in weight per cent and the other oxides in molecular propor- 
tion. The value 47 used comes from determination of the values "a” and ,f b n 
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Fig. 5. Graph showing the spatial variation of 

the 0~value and the volume of the Quaternary 
volcanic materials across the East Japan 
volcanic be t (after Sugimura 1967). 
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for the relation Si0 o = a + b (Na o 0 + K o 0) / A1_CL, Studies in Japanese 

i. I l l J 

Islands show that the isobaths of the seismic plane (mantle earthquakes) 
are parallel to the volcanic front and also to the boundary lines of 
petrologic provinces. This relation is also found in the Aleutian arc, 
though the petrologic boundary is not as clear (Barazangi and Dorman 1969). 

The volume of the volcanic material decreases with the increase in 
distance from the front (Sugimara and Uyeda 1968). It is found that the 
0 - index at the volcanic front is the highest for any part of a given 
island arc (Fig. 5) « For the Aleutian arc it seems to be generally con- 
stant as follows: 

TABLE I 


N. E. Umnak 

36.9 

S. W. Umnak 

39.0 

Adak 

35.3 

Makushin 

37.5 

Elder Point 

36.1 

Risk a 

39.1 

Little Sitkin 

38. 3 

Trident* 

36.88 


(after Sugimura) ,* calculated from (Forbes et al 1967) 

The volcanic front of an island arc may be used as one of the char- 
acteristics of the primary magma for the island arc. 

Gravity Data 

A local negative gravity anomaly has been mapped in the Katmai area 
by Kienle (1968); -20 ^ -35 mgal anomaly is centered on the flank of Mt. 
Trident (Fig. 43) . Kienle (1969) interpreted this anomaly as due to a 
very shallow mass deficiency near. Mt. Trident, the maximum depth to its 
center being about 3 km. The density contrast between andesitic magmas 
and the host rock at Mt. Trident are such that the mass deficiency cannot 
be. due to an andesitic magma chamber. It is interpreted as a volume of 
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host rock heavily impregnated with low density vesicular andesites and 
therefore the gravimetric ally detected mass deficiency under Mt. Tri- 
dent may be related to the chambers discovered seisraically that underlie 
the area at about 10 km depth (Kubota and Berg 1967). Kubota and Berg 
(1967) also' claim a direct relationship between Mt, Katmai and Mt. Nova- 
rupta, however, the proposed chambers have diameters of no more than 5 
km and a 10 km depth. 
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Complete Bouguer (2.67) gravity map of the lit. Trident-Knife 
Peak area. 
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CHAPTER IX 

METHODS AND TECHNIQUES 

Equipment 

Tl\e proton precession magnetometer measures only the total magnetic 
intensity. It was first used to measure the Earth's field by Packard and 
Varian (1954) . The method depends on the fact that in the presence of a 
magnetic field the spin axes of a small proportion of the protons in li- 
quids such as water, become aligned in the direction of field; if the 
direction of the field changes, the protons will precess about the new 

field direction with a frequency given by 211 f - r F where r is the gyro- 

P P 

magnetic ratio of the proton, and F is the field strength in gauss. The 

value of r adopted by the IUGG is r = (2,67513 ± 0.00002) x 10^/0e. sec. 
P P 

The working relation for F is the F = 23,487.4f gammas, (F = 234,874 

—6 

x 10 f gauss). For the Earth's field f lies between 1.5 and 3.0 kc/sec. 

The equipment supplied with a magnetometer includes the sensing ele- 
ment, magnetometer console, interconnecting cable, and their associated 
power supplies. The sensing element is a toroidal coil surrounding a 
hydrocarbon fluid with selected characteristics. This is mounted in a 
container suitable for the proposed application, i.e., a “bird" or a 
"fish” for airborne or seaborne work. The hydrocarbon used in this case 
was kerosene. By energising the coil surrounding the hydrocarbon, the 
protons were aligned in the direction of the applied field. After allow- 
ing a few seconds for the alignment of the protons to reach equilibrium, 
the polarizing field is suddenly ' removed. The aligned protons then pre- 
cess round the Earth's . field , causing an alternating voltage, initially 
of the order of 1 pV, to be induced in the. coil. Tills voltage, which 
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decays as precessional coherence is lost, is applied to a frequency- 
measuring circuit through a tuned high gain amplifier. The minimum 
operating cycle for reliable results is about .3 sec. 

In practice, the coil is arranged co-axially with respect to the 
head, so that it makes a satisfactory angle x^ith the geomagnetic field 
everywhere except on headings close to the magnetic meridian near the 
magnetic equator. 

A system employing two mutually perpendicular coils may be used any 
where. < 

The general arrangement of a magnetometer using this method of mea- 
surement is shoxm in Fig. 8. 

The console operates directly from a source of DC power ranging from 
22 volts to 32 volts. 

The time sequence is shown in Fig. 9, 

One of the portable proton precession magnetometers was used as a 
base station to record the daily variations. The sensor x^as set on the 
top of an 2.4 m foot pole near the base camp. The principle of operation is 
the same as for the above equipment. 

The recorder has a somex^hat small scale and in reading the variations 
the reading errors are about ± 5 gammas. The daily variations at the base 
camp are shox^ in Fig. 13. 

The second portable magnetometer was used for ground surveys, and 
is useful for. accurate investigations of small areas to look at shalloxtf 
structures. 

The airborne survey was carried out using a Hiller 12E helicopter. 

The bird was lifted' from the' ground on the' end of its 80.5 m cable by tak- 
ing off vertically. Landing was done using the reverse procedure . The 
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Fig. 9. Time Chart. 
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Hiller 12E helicopter is small and only allows two people together with 
the magnetometer. 

Methods of Investigation 

Airborne surveys are usually made along a series of parallel tra- 
verses oriented as closely as possible to right angles to the trend of 
the geologic units. The direction may be chosen parallel to some obvious 
geographic feature to aid navigation. The spacing between traverses is 
usually one chosen to be approximately equal to or less than the depth 
of burial of the magnetic source. Perhaps the most difficult problem in 
an airborne survey is to determine the location of aircraft at any instant 
in time. The position of a number of distinct points such as mountain 
tops, cliff edges, and creek or stream junctions were located on maps, 
and in some cases large red panels were laid out for landmarks . 

The landmarks could be clearly seen from the start and end positions 
of each flight line, a direct course between them being flown by eye. The 
beginning and end points, and the crossing of intermediate landmarks were 
also judged by eye, each being marked on the magnetometer recording using 
a hand controlled event marker. Because it is possible to fly slower, sur- 
veys using helicopters are easier than those using fixed wing aircraft, 
however, altitude, speed and direction are not as stable. 

To process the data the magnetometer chart was divided into one- 
tenth mile sections by extrapolation of the known distances between navi- 
gation points on the map, and the length of record between successive 
marks representing these landmarks. The precision of this method is not 
high, and differences can be seen in the same line flown twice at differ- 
ent altitudes, however, as each line was flown at least twice in opposite 
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directions (Fig. 16), and the anomalies are generally very large, this 
is not considered a serious hazard. 

Interpretation of the magnetic profiles depends on ground elevation, 
the flight altitude was recorded manually and the ground elevation 
read off the map for each of the profiles. The flight directions were 
often at different azimuths which were also recorded. Many of the flight 
lines concentrated around the center of Novarupta. 

The starting point and end points were either marked with flags on 
the ground or consisted of conspicuous topographic features, which were 
located using a geodimeter and a theodolite. 

The plotting was done on a 1:21,120 scale map. Because of the errors 
inherent in the method, the cross points of different lines are not well 
defined, however, they are very close. 

Data Analysis 

The total intensity of the magnetic field at any point on the surface 
of the earth can be obtained mathematically by means of a harmonic analy- 
sis of values obtained by direct observation; potential, 
w — a E £ (g m cosm + h™ sinm <J>) P m (cos $) , 
where the values of a, g, h are known, and assuming no crustal sources 
contribute sensibly to the harmonics of degree less than n. 

The magnetic field at the Earth’s surface (Fig. 10) can be repre- 
sented approximately by the field of an eccentric dipole. The difference 
between the observed field and the field generated by a dipole is con- 
sidered as the anomalous or non-dipole field. The anomalies that are gene- 
rated in magnetic materials in the Earth’s crust ai'e called regional ano- 
malies and are assumed to be due to magnetic minerals. These can only ex- 




Fi'j - 10 . The total intensity of the earth's magnetic field expressed in 
kilo gammas . 
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ist up Co the Curie temperatures of the minerals causing them, so 
must be located relatively near the surface of the crust. 

Regional intensity anomalies are often related to magnetite deposits 
or rocks rich in magnetite, and the wave length of the anomaly give an 
indication of its depth. Long wave anomalies usually indicate that the 
magnetic body is situated at greater depths. 

The earth's magnetic field increases in intensity from south to north 
at a uniform rate of 4-5 gammas per mile (Brosge and others 1970). The 
Katmai area has such high local gradients and is so small that it is 
reasonable to neglect this gradient and use the mean values derived from 
the average of all the data as the base value (53655 gammas) . The true 
value should be obtained from a spherical function of the earth's magnetic 
field, but as it is only the relative intensities related to local struc- 
tures that are of interest, the total magnetic intensities used were ob- 
tained from: 

F , - F- - F where F , - observed value 

ob m an ob 

F = mean value 
m 

F - magnetic anomaly 

Any anomalies found indicate a source underground, such as high magnetic 
susceptibility rocks, or non-magnetic materials, the later appearing as a 
relative decrease in magnetic intensity. If we know the regional field, it 
is possible to calculate the magnetic dipole moment of a source capable 
of producing the observed anomaly. In the case of an aeromagnetic survey, 
that calculation is not important, because the magnetic intensity is in- 
versely proportional to the cube of the distance from the magnetic body 
The surveys flown at a high altitude in Katmai give intensities that are 
almost the same as the Earth's magnetic dipole field, except in the 
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vicinity of Novarupta. 

The instrument drift was checked, and daily variation eliminated by 
means of measurements made at the base station. The aeromagnetic data 
were corrected for daily variations, when magnetic disturbances exceeded 
the base value of 53700 ± 5 gammas. 

In the auroral zone magnetic disturbances show large variations with- 
in short distances, though the distances are still large compared with the 
scale of the Valley, so any important variation would have been seen 
at the base station. 

The reduced total intensity anomaly maps shown below were compiled 
from this data using an IBM 360 computer and the program given in Appen- 
dix II . 

Interpretation of Model Studies 

It is possible to estimate the volume and depth of a magnetic body 
from the anomalous field produced at the surface. This will be due to the 
total intensity of the Earth's magnetic field plus that due to the buried 
magnetic material. The intensity of a given magnetic body can be calcula- 
ted as follows (Fig. 11) : 
let F - total intensity observed 

H * total intensity of the regional field 
B « field due to the buried mass 

Using a coordinate system with its origin at the dipole representing the 
buried mass, and letting 

0 - the angle to the observer relative to the dipole axis 
ra - the effective dipole moment 

r « the distance from the dipole center to the observer 
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Fig, 11. Idealized model profile and vector composition 
(see Model Studies p. ). 
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then 


H » h! + H t 
6 8 r r 

B = fit + H f 
8 0 r r 


where 


B = m sin0, 

A 


B - 2m cos Q 

r _____ 


Hg = -Hsin ( 7i — 0 ) 


H = Hcos (it-0) 
r 


F - ( 2m cos0 + Hcos 8)1 + ( m sinQ - HsinG)l 

3 r 3 

r r 


6 


then the total intensity 

F = /(Fr^ + F„^) = /{cos^0( 2m + H)^ + sin^8( m - H) } 

9 3 3 

r r 




where 

h ~ r sin (0-cL) 
therefore 

F = /{ cos 2 8( 2 tr sin^(6 - a) + H)^ + sin^6( m sin^(O-a) - H)^} 


We can get the model profiles from above equation for any given volume, 
magnetic moment and distance. 


Me as urement s 

A series of flight lines was proposed on the basis of the ground 
survey made in 1969 by M. Trible (1971). In practice the lines flown de- 
pended on the weather, and it was not possible to make a systematic sur- 
vey. The aeromagne tic survey was carried out between July 8 and July 28. 

As a result, the flight altitude direction and helicopter speed were 
variable. The flight lines were concentrated in the vicinity of Novarupta, 
but a number of Valley crossings were also made. 
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Fig. 12. Total magnetic intensity profiles near Base Camp (BMC).' 
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Fig . 13. Daily variation of total magnetic intensity on the Baked mountain camp 
(4 feet pole). Base line for each day 53700 gammas. 



A base station was established at Baked Mountain camp, the location 
of which is 58°17.5’N, 155°12.0'W and the elevation is 899.3 m. To obtain 
a value for the magnetic intensity at the base camp, two lines 66 m long 
and 0.6 m apart were marked with stakes at 0.6 m intervals. Magnetometer 
readings were made at each stake, the variation of intensity was 100 
gamma over this area. Measurements were made at two probe heights, 2,4 m 
and 1.2 m from the surface (Fig. 12). The upper values show small diffe- 
rences from the lower ones, the mean values being 53680 and 53700 gammas. 

The aeromagnetic mean value of the base field used was calculated 
from all the available data. Using the average value obtained for each 
different altitude, the individual mean fields were calculated. Fig. 15 
shows the values obtained, and also shows a pronounced difference com- 
pared with total field values obtained from spherical harmonic analysis. 
The mean of the mean values was used as the regional field. The anomaly 
value AF was calculated from: AF = F observed - F mean field — F daily 
variation. The mean field used to prepare the anomaly maps was 53655 y. 
This method is very convenient for looking at anomaly patterns and tec- 
tonic structures over small areas such as the Valley of Ten Thousand 
Smokes. Table 2 presents the direction, altitude, speed, time and course 
of the individual profiles. The isopach maps of geomagnetic field anomaly 
for each group of flight lines indicates the regional residual intensity. 
Six flight levels were used, 550-670ra, 730-760m, 915-975m, 1005“1070m, 
1160- 1190m and 1220m above sea level. The flight altitudes chosen depended 
on the topography and the weather. 

There were 62 aeromagnetic lines flown, of which two- thirds were 
concentrated in the vicinity of Novarupta, This area is important in 
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TABLE 2 


Flight 

Flight 

Elavation 

Time 


Flight 

Flight 

Elavation 

Time 


Course 

Direction 

m, feet 

(Date 

Time) 

Course 

Direction 

m, feet 

(Date 

Time) 

Number 






Number 






2 

SW 

193.5 

745 m(4000) 

July 13 


34 

NE 

37.0 

640 m(2100) 

July 18 

19:50 

3 

SW 

200.5 

760 m(2500) 



35 

SW 

227.5 

760 m(2500) 

July 25 

13:30 

4 

E 

66.0 

775 m(2550) ... 


^ . 

36 

NE 

211.5 

760 m(2500) 



5 

E 

50.5 

760 m(2500) 



37 

SW 

27.3 

1070 m(3500) 



6 

NW 

278.0 

915 m(3000) 



38 

NE 

194.5 

1160 m(3800) 



7 

SW 

219.0 

760 m(2500) 



' 39 

W 

241.5 

1160 m(3800) 



8 

E 

74.5 

610 m(2000) 



40 

NW 

290.0 

1160 m(3800) 



9 

Sw/E 

211.5 

915 m(3000) 



41 

N 

136.5 

1160 m(3800) 


14:30 

11 

E 

50.5 

irregular 



42 

' N 

139.9 

760 m(2500) 


10:00 

12 

W 

50.5 

975 tn(3200) 

July 15 

12:20 • ' 

43 

N 

248.5 

760 m(2500) 



13 

E 

50.5 

1080 m(3550) 



44 

E 

40.5 

775 m ( 2 55 0) 



14 

N 

309.0 

1065 m(3500) 



45 

S 

139.9 

915 m(3000) 



15 

E,W 

254.0 

610 m(2000) 


16:30 

46 

N 


915 m(3000) 






670 m(2200) 



47 

E 

92.3 

915 m(3000) 



16 

S 

153.5 

670 m(2200) 



48 

,W 

250.0 

975 m(3200) 



17 

E 

65.1 

670 m(2200) 



49 . 

E 

75.5 

975 m(3200) 



18 

S 

158.5 

745 m(2450) 



50 

S 


1005 m(3300) 



19 

S 

170.0 

915 m(3000) 



51' 

S 

143.0 

1005 m(3300) 



20 

NE 

23.5 

915 m (3000) 



52 

W 

245.5 

825 m(2700) 



21 

W 

227.8 

915 m(3000) 



- 53 

E 

251.0 

730 m(2400) 


10:50 

22 

SE 

117.5 

915 m(3000). 



54 

N . 

328.5 

790 m(2600) 



23 

E 

72.2 

915 m(3000) 



55 

S 

166.0 

1220 m(4000) 

July 4 

19:30 

24 

N 

345.5 

915 m(3000) 



56 

N 

344.7 

1220 m(4000) 



25 

SW 

194.5 

1190 m(3900) 



57 

S 

17 2.0 

1220 m(4000) 



26 

N 

317.0 

1190 m(3900) 



58 

N 

353.5 

1220 m(4000) 



27 

NW ' 

300.0 

760 m(2500) " 

July 17 

19:30 

59 

S 

169.5 

1220 m(4000) 



31 

NE , SW 35.0 

550 m(1800) 

July 18 

19:50 

60 

E 

79.2 

1220 m(4000) 

July 5 

11:30 




670 m(2200) 



61 

W 

281.5 

1220 m(4000) 



32 

SW 

43.5 

580 m(1900) 



62 

. N 

346.5 

1220 m(4000) 






700 m(2300) 









33 

NW 

49.0 

610 m (2000) 











219.0 

730 m(2400) 
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Fig. 15. Average values of Magnetic total intensity for different flight levels. 
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terms of investigating the role of Novarupta in the 1912 eruption, however, 
preliminary and analysis indicates the need for more data in the surround- 
ing areas . 

The uneven distribution of flight lines makes computation difficult. 
Each of the profiles, Nos. 11, 12, 13, 15, 25, 31, 32, 33, 38, 39, 45, 

51, were flown at two or more levels in both directions. Generally the 
reversed profiles correspond very well with the original, however, as the 
speed was usually uncertain and different in each direction due to wind 
changes, the extrapolation between known markers on the profiles has an 
error due to the speed/distance estimates. An example of the general 
correspondence is seen in profile No. 33 (Fig. 16). The location of the 
profiles are shown in Figs. 32, 35, 38, 41. The individual profiles are 
described below: 

No. 15 E and W (Figs. 18, 29, 31) 

This line traverses the Lethe Valley at 660 m and 726 m. The profiles 
are generally similar, the higher one gives lower values and a smoother 
curve. There is a difference of about 80 gammas between the lines, which 
is comparatively large compared with the other profiles. This may be an 
indication that this area has a thicker flow unit beneath it than those 
further down valley. 

No. 31 (Figs. 19, 29) 

Line No. 31 was flown at 594 m and 726 m. The residual differences 
between the two lines are less than 60 gammas on the west side of the 
valley, reducing so that the whole east side showed no difference between 
the two elevations. 

No. 32 (Figs. 20a, 29) 

This profile has the same pattern as line 33, but shows more fluctua- 
tion in the readings obtained. 


35 







E(BKM) 


-53800 




AS 
y I 


I 5 W 


• I00GAMMAS 


6!0M 


670M 


0.5 KM 


- 540 m 



Fig. 18. Total magnetic intensity profiles for two flight levels 
in the Lethe Valley (see Fig. 25, 33). 
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Fig. 19. Total magnetic intensity profiles for two flight levels in the 
Main Valley (see Fig. 29). 





No. 33 (Figs. 20b, 33) 


This line was flown at 792 m and 726. ra, and had an upper limit of 
the be tween-line differences of about 50 gammas, although there was an 
elevation difference of more than 66 m. 

Lines Nos. 12, 13, 39, 11 and Ground Traverse (Fig . 17) 

These lines were centered around Novarupta, and the intensity decreas- 
ed with increasing elevation of the flight lines, the vertical gradients 
also decreased, being about 1,14 gammas /m between 1160 m and 1065 m 1.33 
gammas /m from 975 m to 915 m. The reason for this may be explained by the 
higher magnetisation of the dacitic rocks forming Novarupta Dome. Thin 
section studies on these rocks show about 3% magnetite. 

The ground survey profile (Fig, 14) in the same area showed many 
more small scale variations, the highest of which was a 1450 y anomaly near 
the summit of Novarupta. All the profiles showed a positive anomaly lo- 
cated around the lava dome, which is presumably due to the remnant magne- 
tization of the rocks composing Novarupta Dome. 

Lines Nos. 41, 45, 58, 59 (Figs. 24, 25) 

These lines are located on Katmai Pass and on the flanks of Mt. Tri- 
dent, crossing Novarupta and Trident Basin to Broken Mountain. Parts of 
these profiles shown an inverse relationship between relative land elevation 
and magnetic intensity, which implies that the lower flanks of Trident 
and Broken Mountain are composed of low susceptibility Naknek sediments. The 
magnetic difference is clearly seen in the profiles. The vertical 
gradient is + 3 ^ .8 garamas/m. 

Line Nos » 25 ' and 38 (Fig. 21, 39) 

'These profiles traverse Novarupta from south to north at 1155 m and 
1187 m. They show symmetrical magnetic features all positive with respect 
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Fie. 20 b. Total magnetic intensity profiles at two flight levels across tne 
Main Valley (see Fig. 29, 30). 





to the base line of 53655 y. 

The anomaly pattern is about three times as wide as the visible 
part of Novarupta lava dome. 

No. 56 N (Figs. 22, 42) 

56 N (1220 m) crossed from Falling Mountain via Novarupta basin to 
Baked Mountain. The part of the profile across Falling Mountain is typi- 
cal of an andesitic volcanic feature. The maximum intensity is 56920 gam- 
mas. The steepness of the gradient at the base of Baked Mountain indicates 
a fault structure between the Naknek sediments of Baken Mountain showing 
low intensities, and presumed volcanics in the Novarupta Valley. 

No. 57 S (1220 m ) (Figs . 23, 42) 

This profile crossed Broken Mountain via Novarupta to Katmai Pass at 
1220 m. It is typical for a body magnetized in the direction of the present 
field. The depth of an equivalent magnetic dipole has been calculated 
(see p ). Stumbling Hill is of smaller intensity but similar to Nova- 
rupta. The relation of Stumbling Hill to the sediments of Broken Mountain 
is not clear from the profiles. 

No. 61 W (1220 m) (Fig. 42) 

This section crossed Trident basin via Novarupta, Novarupta basin, 
and the Lethe Valley. Compared with the lower altitude flights, the sec- 
tion over Novarupta gave magnetic intensities of low frequency and low 
amplitude. 

The section crossing the main valley looks almost flat indicating 
that this pyroclastic flow has a low effective magnetic intensity. 

No. 62 N (Fig. 42) 

This flight line, from the base of Mt, Mageilc along the west flank 
of Baked Mountain and across the Knife Creek branch at 1320 m gives an 
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Fig. 23 . Total intensity aero-magnetic profile over Novarupta. 
Flight level 1220 m (see Fig. 42 ). 
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almost flat magnetic profile and shows no indication of underground 
structures or magnetic bodies . 

Ground Traverses (Fig. 14) 

Although the aeromagne tic anomalies do not give any indication that 
the near surface material is irregularly magnetized, a very detailed 
magnetometer survey of the outcrops revealed a highly irregular small 
scale anomaly pattern. This indicates that the magnetization is not uni- 
form over distances of more than a few meters. 

The ground magnetic survey covered the head of the Valley of Ten 
Thousand Smokes (Figs. 26, 27, 28). Two of the survey lines corresponded 
exactly with seismic profiles, these being along a median line down the 
Intermediate Valley, between Baked and Broken Mountains (Fig. 27), and a 
median line between Falling and Baked Mountains in Novarupta Valley 
(Fig. 26). A ground traverse was also made across Novarupta Valley on a 
line through the peaks of Falling and Baked Mountains (Fig. 17) * The lo- 
cation of the magnetometer stations was obtained using a tape measure; 
these data are corrected for level changes, i.e., corrected to the hori- 
zontal distances between stations. The elevation of each position was 
measured using a barometric altimeter. As the magnetometer used for these 
surveys normally served as the base station, no record of the daily varia- 
tion was available. The readings at the base station before and after the 
surveys were noted, and the change over the few hours involved, was assum- 
ed to be linear. The data were plotted on 1:21120 scale maps. The length 
of sensor support pole was 8 ft. for all these surveys and care was taken 
to have no magnetic material around the probe. 

The most significant features are associated with fumaroles, this 
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was particularly true for the line along the axis of Novarupta Valley, 
where a line of fumarolic holes are located near the foot of Nova- 
rupta. The fumaroles are not currently active. The relationship between 
the fumaroles and the high magnetic anomalies is not clear, but is pro- 
bably associated with the deposition of magnetite. Near the foot of Nova- 
rupta a depression about 20 or 30 m across contains many volcanic bombs 
and volcanic breccia. It is possible that this ejecta might also be re- 
lated to the observed magnetic intensity. These profiles showed generally 
high amplitude and short wavelength anomalies (see also Tigs. 17,26), 
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CHAPTER III 


MAGNETIC DATA 

The contour maps used in the following sections were produced using 
the IBM 360 computer. The aeromagnetic data from each range of flight 
levels has been computed separately, and is presented in three maps, an 
index map, a total intensity map and an anomaly map using 53655 gammas 
as a base line. 

In general, the two maps have been plotted with different contour 
intervals, so show features with different scales. It should perhaps be. 
noted here that the computer tends to n ere ate anomalies" for areas in 
which there is no data, so great caution must be used in interpreting 
these maps. 

Since the base line from which the anomaly maps are drawn is the mean 
value of the field within the. area of the Valley of Ten Thousand Smokes, 
the anomaly maps are particularly useful in looking for variations with- 
in the ash flow. 

Novarupta Area 

The rim of the Novarupta structures appear as an oval feature elon- 
gated east— west with the dome of Novarupta towards the western end. 

Xf Novarupta was the sole source of pyroclastics , then the vent of 
Novarupta must have been higher during the eruption to account for the 
present pyroclastic distribution. In this way one eruptive center could 
account for the whole pyroclastic flow. 

The distribution of magnetic anomalies obtained from the 915 m and 
975 m flight altitudes (Figs. 36, 37, 38) showed an oval anomaly pattern 


51 



In the vicinity of Novarupta and Stumbling Hill, which perhaps indicates 
a simple structure. 

The oval magnetic anomaly around Novarupta is not characteristic of 
crater faulting, because the horizontal magnetic gradient is too shallow 
(Fig. 44). A negative magnetic anomaly is associated with the crater rim 
which is probably due to underlying Naknek sediments. In contrast, Fall- 
ing Mountain has NE-SW fault evidenced by the 100-150 m high north face. 
The magnetic profiles over these faults show very steep gradients, and 
Falling Mountain has a very high anomalous field indicating that this 
mountain is of different material, and has a different susceptibility from 
that surrounding it. 

The postulated collapse of ancestral Novarupta probably took place 
shortly after the main eruption, and the position and shape of the elon- 
gate crater was controlled by the overall faulting in the volcanic 
area . 

The changes in intensity with increasing altitude over Novarupta 
dome indicates that the magnetic source is relatively shallow. 

The magnetic map at the 915 m - 975 m flight altitude over the top 
of Novarupta shov/s a conspicious double peak (profiles 5, 9, 12, 3.3) (Fig. 
17) while those flown at 1160 v 1190 ro show a smaller difference of mag- 
netic intensity (Fig. 21). The magnetic contour map gives a pattern simi- 
lar to that expected if the magnetization were in the same direction as 
the present field. Profile number 11 (Fig. 17) crossing Novarupta and 
Stumbling Hill seems similar to that predicted by Rikitake and Hagiwara 
(1965) for a double cone mode (see appendix). They show an anomaly over 
closely spaced cones, from which it is suggested that the two positive 
peaks were probably made by adjacent intruded or extruded material. In 
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profile 11, the second peak is associated with the many collapse features 
of Stumbling Hill, and may thus point to some kind of double volcano. The 
magnetic map based on 1220 m profiles shows only one peak, the double 
peak having disappeared. This indicates that the surface of the magnetic 
body may be roughly concave between Novarupta and Stumbling Hill. In this 
way lower altitude flights would see a double cone similar to those of 
Rikitake’s model (see appendix II), however in high altitude flight maps 
the two peaks in the magnetic distribution can be neglected. The center 
of the second magnetic peak is displaced towards the northeast from the 

center of Novarupta dome (Fig. 40) . 

The concealed magnetic body causing the main anomaly may be the same 
material as the plug of Novarupta, and may be still moving beneath Stumb- 
ling Hill as evidenced by very fresh looking fault scarps. There is also 
the possibility that Stumbling Hill may be a dome similar to Novarupta. 
Curtis (1968) demonstrated subsidence structures which took place after 
the eruption of Novarupta, and estimated an elevation change of approxi- 
mately 150 m. This subsidence alone would not be enough, but the possi- 
bility of removal of part of ancestral Novarupta by explosive activity 
exists . 

The Novarupta phenomena seems to be similar in character to Usu Dake 
in Japan. Nuee ardente type eruptions took place on Usu Dake in 1822, and 
again in 1910, when the Usu appeared as a crypto-dome (volcano) (rising 
daeite magma?), and in 1944-45 it made a lava dome as pseudo-belonite . 

These rocks were dacites and had a high viscosity of the order of 
10^' poise, and a temperature of 980°C (Minato et al. 1965). 

The only large-scale variations of total intensity seen in the aero- 
magnetic surveys were those associated with Novarupta dome, falling Moun- 
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tain and Cerberus. The investigation has proved useful for looking at 
the magnetic effects of the central vent, but surprisingly shows nothing 
that can be clearly interpreted as magnetic peaks away from, but still 
directly associated with, the dome. For comparison, one can consider the 
case of Hakone caldera in Japan, where good agreement is found between the 
magnetic effects and the position of both the cone and the caldera rim. 
This difference may be due to the number of eruptions of the volcano. Ha— 
kone has a more complex history and represents the accumulation of many 
eruptions. The triple volcano of Hakone is 10 km in diameter and later 
activities have been controlled by the first eruption which produced local 
faulting. 

The dominant magnetic trend in the Katmai area does not parallel the 
lines of the Earth's main field, but seems to be affected only by the 
overall geologic structures and hence only indirectly by previous volcanic 
activity. 

The peak of the anomaly pattern seems to be a little to the east of 
the topographic high. In the ae romagne tic profiles, Novarupta shows a re- 
markable positive pattern of about 800 gammas over the dome, furthermore, 
negative anomalies of about —250 gammas appear in a symmetrical pattern 
about a ENE-WSW axis, which is parallel to the tectonic trend. These ne- 
gative anomalies are assumed to be due to sedimentary rocks of low magne- 
tic content (Naknek) . Comparing this map with the 1160 v 1190 m flight 
altitude map (Fig. 41), it seems that the positive anomalies have de- 
creased and are smoother at higher elevations (Fig. 44). The almost cir- 
cular negative anomaly around Novarupta dome derived from the 915 'v 9 75 
m lines (Fig. 38), shows a distribution which also agrees with the sur- 
face geology. The main positive anomaly reflects the exposed lava dome. 
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The. total intensity map derived from profiles flown at 1220 m (Fig. 
42, 43, 44) show features characteristic of the main tectonic structures 
in a northeast to southwest direction. These are offset in a small area 
along an east-west axis making an en echelon structure. Two steep anoma- 
lies are associated with Falling and Cerberus mountains (1165 m, 1080 m) . 

Model Profiles 

From the equation used for the model studies (see Model studies), 

F = /{cos^6( 2m sin^(9-a) + H) ^ + sin^0( m sin^(6-a) - H)^} 



Then, the dipole moment can be represented by: 

m = M V, where M is the magnetic moment per unit volume and V is the 
o o 

volume of the body. 

h = the distance from the dipole center to the flight line, 
a = the inclination of the ambient field. 

Assuming this model, the field expected for an ellipsoidal tube was 
calculated. The values chosen were 

-3 

x - 720m, y = 650m and M - 2 x 10 emu/cc 

o 

where x and y are the semi -major and semi -minor axes of the ellipse. These 
dimensions are based on the observed size of Novarupta and Stumbling Hill 
(Fig. 21 ) . 

In order to obtain profiles using the above relations, it is neces- 
sary to choose values for the three variables V, m and h. This has been 
done using the "best guess" values based on observations. 

From the top of the magnetic body to the flight level was 300m, so 

that 

h = 300m + L/2, where L ® the length of the tube. 
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V = 7Txy(300 + L/2) 
m = 7Txy(300 +L/2)M 

For any given value of h, it is possible to draw a profile corre- 
sponding to the angle 8, The result of these calculations are shown as 
profile in Fig. 32 . The assumptions used are very rough, however, com- 
paring the observed and calculated data it is possible to recognize simi- 
lar features in both. 

If the dimensions of the system chosen are approximately correct, 
then the dipole center appears to be located about 1 km below the 1.22 km 
flight line, or about .9 km below the surface. 

Using the simple model described above it is not possible to make 
both the values of maxima and minima coincide, and at the same time have 
similar values for the horizontal distance between them when comparing 
observed and calculated profiles. A better fit could be obtained by using 
a more complicated shape for the body, and/or increasing the inclination 
of the magnetization. However, the general fit of the two profiles is 
sufficiently good to indicate a relatively shallow source at about 900m 
below the surface for the Novarupta-S tumbling Hill anomalies. 

The exposed lava dome is about 300 m in diameter, so if the assump- 
tions are correct, then there is more magnetic material on the northeast 
side of the present lava dome. 

A similar calculation using the technique of Vacquier et al (1951) 
gives a value for the long axis of the magnetic body beneath Novarupta as 
about two-thirds of a mile and about half a mile for the short axis. 

Profiles 57S and 61E taken in north -south and west-east direction 
respectively (Figs. 17,23) crossing Novarupta Dome show very elliptical 
magnetic body features. They show large positive anomalies that coincide 
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Fig. 29. Flight lines of elevation 550 - 670 in above sea level. 















Fig. 32. A comparison of various calculated magnetic profiles with the observed 
profile over Novarupta. 




with the rhyodacite lava dome. The dome is composed of highly magnetic 
material but must have a wider distribution than is indicated by the ex- 
posed dome itself. The maximum anomaly on Falling Mountain is 2500 gammas, 
it has a very steep gradient and shows a circular cone-like feature. Using 
similar calculations to those used for Novarupta, the magnetic source has 
a diameter of about 700 m and the dipole center is located at 1000 v 1050 
m below the flight line. These high gradients for the magnetic intensity 
and the shallow dipole source, support the idea that the volcanoes of the 
Katmai area are younger than Tertiary. This is because the data precludes 
a thick sequence of volcanics, unless the temperatures beneath the vol- 
canoes are higher than the Curie Point, also the magnetic distribution 
would be more complex. 

Intermediate Va lley and Fumaroles 

Magnetic anomalies of fumarolic origin were not generally seen in the 
aeromagnetic. maps except for those in the southern valley located on the 
supposed fault line and along the flank of Baked Mountain (Fig. 35). This 
area of fumarolic activity appears to be relatively large and Involve a 
large volume of magnetic material. It is possible that the acidic gases 
and heat might be supplied by some kind of intrusive source or magma 
chamber. In general, there appear to be two kinds of fumaroles: a) those 
that originate in the pyroclastic flow (nuee ardente) , and b) those of 
deeper origin from at least underneath the Naknek sediments which are 
fault controlled (Gedney 1969) . 

Ground magnetic profiles along the center line of the Novarupta Val- 
ley (Fig. 17) show magnetic anomalies that correspond exactly with aero- 
magnetic profiles, except for a few sharp, and presumably very local, ano- 
malies. The narrow magnetic anomalies of high frequency and high amplitude 
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in the vicinity of furaaroles are able to be directly related to the 
surface expression of fumarole activity (Trible 1971). The high ano- 
malies also indicate a source which is both shallow and narrow. 

The ground magnetic profile which traversed the Novarupta Valley 
from Falling Mountain to Baked Mountain also showed fumarole associated 
anomalies (Fig. 25) . 

Intermediate Valley 

The Intermediate Valley has a thin pyroclastic flow in comparison 
with the other branches of the Valley (Kienle 1967, Fenner 1922) . A de- 
tailed survey over a fumarole area (Trible 1970) again gave high frequen- 
cy and high amplitude anomalies, 

A profile along the median line of- the Intermediate Valley was ob- 
served at two heights, 2.4 m and 1.2 m (Fig. 27). If the data from these 
two elevations are plotted, it can be seen that the sharp high amplitude 
anomalies also have high vertical gradients of total intensity of about 
16.7 ^ 2,5 gammas /m again indicating a shallow origin. The higher intensity 
anomaly areas are all normally magnetized, but the magnetization of the 
other smaller anomalies are uncertain. 

The margin of Intermediate Valley . shows one of the highest anomalies 
(250-300 gammas) recognized in the 730-760 m aeromagnetic survey. This 
magnetic high correlates x^ith the fumarole grid of Trible (19 70) and in- 
dicates a considerable mass of magnetic material. According to Trible 
(1971), even where large furaaroles existed, the magnetic minerals concen- 
trated by the fumarole would leach and weather out very quickly after acti- 
vity stopped, Trible also suggests that it is possible for the magnetic 
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minerals to concentrate in the mouth of the fumaroles and to occasionally 
be preserved. Where the magnetic minerals are preserved, a magnetic ano- 
maly exists, but the absence of magnetic expression does not preclude an 
old hot fumarole. Therefore, it is suggested that the fumaroles along the 
margin of the Intermediate Valley were active until comparatively recent 
time, and the magnetic minerals were formed hydro thermally in the mouth 
of the fumarole vents, and were also concentrated at depth and presumably 
preserved in some way. 

Knife Creek, Lethe River and Main Valley Profiles 

Seven flight lines in the Main Valley and Lethe Valley were flown at 
550-670 m and the results contoured (Fig. 29, 30, 31). The total field 
measurements showed flat profile when compared with those from the head 
of the valley, but with both sides of the main valley and the northside 
of Baked Mountain showing lower intensities. In other words, the low 
values appear to correspond with the sedimentary rocks (Naknek formation) 
which contain only a small amount of magnetic material. 

The anomalies were mostly less than 250 gammas with respect to the 
average value. This small value could be caused by a source within the 
pyroclastic rocks. 

The contour map of the 730-760 m flight line data (Fig. 33, 34, 35) 
shows four positive and three negative anomalies in the upper part of the 
valley. The down valley portion showing the same situation as the 550-670 
maps, but with the values decreasing in accordance with a vertical 
gradient of .2 gammas . The anomalies over the valley floor are not larger 
and are probably due to the irregular topography of the U-shape valley 
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left by the last glacial retreat. Significant anomaly patterns seen by 
contouring the data are negative and largely associated with areas of 
Naknek sediments. A positive anomaly extends over a comparatively large 
area in the Knife Creek Valley, a portion of Broken Mountain and the 
lower end of the Intermediate Valley, however, in the center of the posi- 
tive area is a negative anomaly. This negative anomaly may be a function 
of the line spacing, but the overall positive anomaly appears real. The 
steepness of the horizontal gradients point to a fault structure, and may 
well be connected with the higher fumerolic activity of the area mentioned 
before. 

With regard to the Main Valley, the total intensity decreased as 
expected with higher elevation, i.e., the gradient of the magnetic 
field probably represents a normally magnetised feature. The smoothing 
of the curve shows that the sources of the irregular field are shallow, 
and these profiles show nothing of the deeper structures (Fig. 31, 35). 

Profiles in the Intermediate Valley at 1160 and 1190 m, (Fig. 39, 40, 
41) suggests a complicated pattern in the distribution of thickness of 
volcanics, and by implication a complicated paleotopography . 

Negative anomalies correlate with the head of the Intermediate Valley, 
and the Novarupta and Trident basins, however, these anomalies show com- 
paratively lower amplitudes for higher flight levels as might be expected 
for a relatively shallow source. At the 1190 m flight level the negative 
anomaly associated with the Naknek sediments of Baked and Broken Mountains 
is reduced. This is presumably because the magnetic bodies causing the 
positive-negative contrast are shallow and thin. At the higher altitudes 
the magnetic profiles do not show any marked effect directly attributable 
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Fig- 40. Total intensity aero-magnetic map over Novarupta, flight height at 1160 
1190m. BRM: Broken Mt.. BKM: Baked Mt. 



















to the pyroclastic valley fill.' This fact suggests that the pyroclastic 
flow is not very thick. 

Pyroclastic Temperatures 

The minimum welding temperature of acidic volcanic glass is about 
600°C (Boyd 1961). It is also possible to obtain an idea of the tempera- 
ture regime from the magnetic properties of the rocks. The presence of a 
finely dispersed uniformly magnetized magnetite must have been formed 
under relatively high temperature condition (above the Curie Point for 
magnetite (560°C)). 

The welding and deformation of glass particles can take place at a 
very low pressure, but the temperature, necessary for a complete welding, 
should be at least 800°C (Smith and others 1958). Lovering extrapolated 
back the measured temperature of fumarole No. 1 measured by Zies at 600°C, 
and found it to be between 800°C and 900°C. The fact that the pyroclastic 
flows were found to be uniformly magnetized (Trible 1971) along the 
direction of the present magnetic dipole field, and were magnetically 
stable, demonstrated that the temperature was higher than the Curie Point 
of the minerals involved (Packer, personal communication). The pyroclastic 
flows sampled were not all welded, which indicates lower temperatures 
prevailed for some areas for some of the time, or alternatively that 
some areas were heated more than others after deposition. 

As shown above, the pyroclastic flows, as opposed to the interbedded 
airfall deposits, must have been laid down at temperatures between 560° 
and 800°C . 

The average susceptibility of the pyroclastics from the Lethe indu- 
- 1 

rated ash is 1.33 x 10 emu/cc (Trible 1970), with a maximum value of 
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1.76 x 10” 3 emu/cc. Trible (1971) also noted that - "examination of the 

character of the tuff comprising the individual discs disclosed the 

correlation of colored streaks with measured susceptibility: white and 

brown streaks are prominent in discs of lower susceptibility, black bands 

are associated with the higher values". The Knife Creek pyroclastic flow 

—3 

units have a mean value of 1.8 x 10 emu/cc and other susceptibility mea- 
surements on the tuff unit from various locations in the Valley give com- 
parable values. 

Proven pre-1912 pyroclastic flows have not been observed and the 
similarity of the susceptibility of the different tuff units would indi- 
cate that if they exist, they had a similar origin to the 1912 flow units. 
Thus, pre-1912 units would not be distinguishable magnetically. 

Samples of the Novarupta banded lava gave a bimodal distribution of 
magnetic susceptibility values. The gray glass gave a susceptibility of 

O 

0.3 x 10 emu/cc in contrast with an average for the prominent dark 
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streaks of 1.4 x 10 emu/cc (see Table 3). 

Trible concluded that both the Naknek sediments and the fumarolic 

altered ash can be represented by 0.5 x 10 emu/cc. The susceptibilities 

of the airf all pumice, the water-laid tuff, the indurated ash and the 

-3 

higher mode Novarupta lava gave 1.2 x 10 emu/cc. Hydro thermally altered 
sediments and pyroclastic rocks often gave high susceptibility, presumably 
due to deposition of magnetic minerals. In many cases, as the fumarolic 
activity becomes weaker, the hydrothermal activity changes in character, 
and resolution of the magnetic minerals takes place. Unless protected by 
an impervious layer, the high intensity magnetic minerals disappear. 

The altered air fall tuff sediments have a high susceptibility in 
some fumarolic areas, as indicated by the magnetic ground surveys, though 
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the volume of effective magnetic minerals is limited. When the elevation 
of the survey lines was higher than the magnetic mass diameters, their 
effect was very small. 


TABLE III 


Mt. Novarupta 

Dacite 

Dark 

band 

-3 

1,4 x 10 emu 

Gray 

G1 ass 

0.3 x 10~ 3 

Naknek Formation 




1.0^6.5x10 ^ 

. 

Indurated Tuff 




1.2^1.3x10' 3 
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CHAPTER IV 


CONCLUSION 

The most outstanding feature of the magnetic anomaly pattern is the 
magnetic high found in the vicinity of Novarupta. This suggests that the 
rocks underlying this area are more strongly magnetized than others In 
the surrounding area. The magnetization of the rocks composing the ex- 
posed dome of Novarupta are generally not very strong compared with other 

-3 

volcanics, having susceptibilities around 1 x 10 emu/cc. 

The magnetic anomaly pattern of Novarupta as seen from lower eleva- 
tions showed high frequency and high amplitude patterns such as would be 
generated by double volcanic cones with shallow structures. However, these 
small scale anomalies disappeared at higher elevations and indicate a 
single lava dome. 

The crater rim features and moat visible around Novarupta are re- 
flected in the magnetic anomaly maps by low amplitude and low frequency 
magnetic patterns, these are typical of those observed over thick Naknek 
deposits. The obvious interpretation is in terms of the surrounding ma- 
terial being low susceptibility sediments, through which the plug or dome 
of Novarupta has been pushed. 

The oval anomaly patterns around Novarupta are elongate in an ENE- 
WSW direction, and include Stumbling Hill. Although no outcrops were found 
on Stumbling Hill, there is good reason to think that this hill is formed 
by intruded or domed— up material similar to that of Novarupta, 

The magnetic anomaly pattern changes suddenly at the edge of Falling 
Mountain, Mt. Mageik and Mt. Trident. This is presumably due to the mass 
of magnetic material forming these mountains. 
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A pair of negative and positive magnetic anomalies is associated 
with each of Novarupta, falling Mountain and Mt. Cerberus, and is similar 
to that expected if each of these is considered as a dipole magnetized in 
the direction of the present geomagnetic field. It is of interest to note 
that the elongation of the Novarupta magnetic anomaly and the volcanoes 
of the area are all roughly aligned parallel to the structural trend of 
the Alaska Peninsula. 

The most acceptable model based solely on magnetic data suggests 
that Novarupta and Stumbling Hill are pipes of magnetic igneous material 
intruded through the essentially non-magnetic Naknek formation. This then 
makes it difficult to visualize a much larger ancestral Novarupta, high 
enough to be the source of the pyroclastic fill for all parts of the 
Valley of Ten Thousand Smokes at once. This would require that either the 
material of ancestral Novarupta was largely low susceptibility and thus 
similar to the Naknek, and the destruction of the original peak was by 
collapse, or, that it was a composite cone similar to the volcanoes seen 
in the area today, and the material was removed explosively. This would 
presumably require that the whole composite cone down to Naknek bedrock 
be removed. Neither of these possibilities fit well with the field evi- 
dence. 

Subsidence by some form of block faulting is also a possibility, but 
these will only be seen magnetically where there is a sufficiently large 
susceptibility contrast. High horizontal magnetic gradients indicative of 
faulting are seen along the base of Falling Mountain, and also associated 
with the northern part of Broken Mountain. The gradients in these two 
areas may well have been enhanced by fumarolic action along the proposed 
fracture zones , 
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Over the flat portion of the Valley, no large marked anomalies are 
found, though a number of small scattered local magnetic highs are usually 
associated with fumarolic features. At the foot of Broken Mountain and 
across the lower end of the Intermediate Valley, negative anomalies of 
the order of a few hundred gammas are found around the recently active 
geothermal areas. 

A number of weak magnetic highs seem to be associated with sharper 
peaks over apparently thicker pyroclastic flows. As these anomalies 
amount to only a hundred gammas or so, and are local in extent, they 
probably represent minor variation in the thickness of flow units or per- 
haps an accumulation of magnetite from fumaroles . 

Trible (19 71) discussed the possibility that the glacial drift has 
a high susceptibility and could account for some of the magnetic highs, 
and it is also possible that some material interpreted as high susceptibi- 
lity glacial drift may in fact be pre-1912 flow units. 
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Appendix I 


Steinberg and Rivosh (1965) concluded that the magnetic anomalies 
were controlled by the shape of the volcanic pile by using the Rikitake 
model. This was especially true for dormant volcanoes, where there is a 
magma reservoir maintaining temperatures over the Curie Point, as it is 
recognized that the anomalies are dependent on the remnant magnetization 
of rocks. Figs. I, II, III and IV show the analytically computed total 
magnetic intensity anomaly at five different elevations above Magnetized 
Circular Cone (after ’ Rikitake and Hagiwara, 1965). 
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Fie. I. AF over a magnetized circular cone, which ^ 
is shown by the shaded circle, at a height of 0.2 in 
units of the bottom radius. The numerals indicate tie 
field intensity when the intensity of magnetization 
and the. bottom radius are taken as unity. It is assume 
that the direction of magnetization lies m the magnetic 
meridian plane and that its inclination agrees with the 
geomagnetic dip which is taken as 48 . The slope ang 
is assumed as 10°. 
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Magnetic Anomaly over a Magnetized Circular Cone 


r 



F 


E 

D 


c - 



Fig. II. AF profiles along the north-south line passing 
through the centre of a cone of which the slope angle is 
10°. The direction of magnetization and its relation to 
the geomagnetic field have been given in the caption of 
Fig. I. The intensity of magnetization and the bottom 
radius are both taken as unity. The flight heights are 
taken as 0.1, 0.2, 0.3, 0.4, 0.5 and 1.0. 
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Jlfapnctic Anomaly over « Magnetized Circular Cone 



Fig. III. AF over two cones closely placed. The line 
connectxng the centres runs in a north-south direction 

iSne“*«i1» S a °d f T COneS 3nd th£ relati - between the 
gnetization and the geomapnetir fi’pirt *.1 

those for Fig. I. ale the same as 
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Appendix II 

nor. fortran iv 360n-fo-679 3-1 prcjfyl2 


DATE 0 3/ 03/71 TIME 01-30.50 PAGE 0001 


C PROGRAM TO DRAW PROFILES FOR MAGNETIC DATA AND ALTITUDES WHEN NEEDED 
C 


. 0001 

0 I MENS ION 

I SYM ( 5 ) XP(.IO.O) 

.YP( IOOJ.t 2 P (.100 ). 



0 o 0 2 

REAL 

NAME! 2) 



0003 

DATA 

I SYM / 0, 12* 5* 

2 *6/ 


0006 ... _ .... . ... 

. C IN = 2.5 7 

.. .... _ _ _ 

... 



c 

C INITIALIZE PLOT TA^E 

0005 CALL PLOTS T 
C 

C READ FIRST CARD OF SERIFS - MIN X,Y, M A X X , Y , NBR PROFILES* HEADING 

~" r C AND RANGE OF ALTITUDE IF NALP IS SET 

C XMIN AND YMIN = STARTING POINT AMO XMAX YM AX = END POINT XMIN MAY BE 

C A 0 r ; j a f |. Y MU ME P. 1 0 ALL Y L A P. G F R THAN XMAX 

C CURRENT LIMIT IS FIVE PROFILES PER PLOT. THERE IS I FILF PER PROFILE 
C. CHANGE r»LN COLORS FOR DIFFER T: N T FLIGHT ALTITUDES IF DESIRED/ 

C ... 

0006 ' 5 iiFAl) 11,1000, F N D = 000 i HEAD, XMIN* YMIN* XMAX, YM AX * 

l AMIN, A MAX , NPR , NALT 

0007 1000 FfPMAT ( 7F10.0, 212.) . 

C 

r ft THAN FIVE (55 PROFILES ARE REQUESTED - TRUNCATE TO FIVE AND 

C SET CC1 )C TO TERMINATE JOB AFTER THOSE FIVE ........ 

C 

QC08 I TERM - 0 

00 00 I F ( N MR _ . LE * 5 1 . SO..TO _ 6 : 

0 0 1 0 M MR - 5 

0011 [TERM - 1 

C CALCULATE LENGTH Hr AXIS, DRAW & ANNOTATE 
C 


0 012 

6 D X = 

XMAX - 

XMIN 


001 3 

QY = 

Y M A X - 

YMIN 


00 16 

XLEN 

= SQR T ( 

D X * OX + 

DY * DY 

GO 1 5 

CA1 1 

PLOT ( 

1.5* 2.0, 3 

) 

00 ! 6 

XL ~ 

1.5 + ( 

XLEN / UN 

) 

0 C 1 7 

CALL 

P! 0 T ( 

XL, 2.0, 2 

J 

0 0 1 ?j 

no \ 

0 I - I , 

s 


00 10 

FAC 

» ( I - 

l. ) / 6. 


GUO 

XU = 

X MIN + 

FAC * OX 


00 71 

YN = 

YMIN + 

FAC * DY 


0 0?? 

X = 

1.5 v FAC * ( XLEN 

/ CIN i 


0023 Y = 2.0 

no PA CAL l. SYMBOL ( X, Y T *16* 13, 0.0, -1 > 

0 0 ? 5 -- Y - l.P 

0 0 ? 6 CALL NUMBER I X, Y, .16, XN, 270.* 2 ) 

0 0^7 CALL SYMBOL (-0.0* - 0 . D , -Q. 0 270. * 

0 0 2 >3 C A L L ' N U M A c R ( — 0.0* —0*0, —0.0* Y N » 2 7 0 * * 2 ) 

0029 10 CONTINUE 

C n R A w A N o L A ■> f L HA G N F: T I C AXIS 




OGS FORTRAN IV 360N-F0-A79 3- L 


PROFYL2 


DATE 03/03/71 


TlMt 


01.30.50 


PAGE 


0030 


CALL PLOT I 1.5, 

0071 


CALL PLOT I 1.5, 

00 32 


YN = 52000. 

00 33 


X = 0.5 

00 3% 


Xl= 1.5 

00 35 


Y = - . 5 

0036 


00 15 I = 1 , A . . 

0037 


YN = YN + 1000. 

00 38 ' 


Y = Y * ? . 5 

00 79 


CALI. PLOT ( X, Y 

00 AO 


CALL NUM7 r R l X, 

COM 


CALI SYNOHL ( XI 

OOA2 


15 CONTINUE 



C 



C DP AW AMD LA3LF ALT IT 



C 

00 A3 


IF I M AIT .EO. 0 

0 0 A A 


OA * A MAX - AMIN 

00 A 5 


ML - DA / 500. *■ 

00A6 


CALL PLOT < XL, 

0 0 A 7 


CALL PLOT ( XL , 

0 0 A 3 


YN = AMI N 

00A9 


X -- XL + .2 

0050 


DO IB I = 1, ML 

0 0 3 1 


F = I “ 1 

005? 


Y N = A M IN + F * 

no 63 


Y = I YN - AMIN 

0 0 5 A 


CM, L PLOT t XL , 

00 55 


CALL S v ? 0 L ( XL 

0 0 56 


CALL MLC'DER [ X, 

0 0 3 ’ 


L R CONTINUE 

0 0 63 


CALL SYM30L { XL 

00 59 


CALL MUMPER { X, 

0080 


If? [ p = 0 



c 



C FND FILE on PLOT TAP 



C 

0 0' 6 1 


? 0 IP - I P + 1 

0 0 6? 


IF < IP . n ,) . 1 ) 

0 0 6 3 


IF < IP .LE. NPF 

00 6 A 


X N F 3 - A MAXI ( > 

C 065 


XNFW = XL FM / C 

0066 


YNFV! = 0.0 

0 0 6 7 


CALL PLOT ( X ‘ ! E r 

0088 


IF ( I TERM , EG. 

0069 


00 TO 5 

0 0 70 


?5 CALL PI.DT [ 0.0 



C 



C R F AD CAPO FOR INOIV 



C 

00 71 


?6 R. F AO (1,10011 AL' 

0 C 7 2 


100 1 FORMAT < FI 0.0, 

0 0 7 3 


T * 0 


0002 


93 ; 


DOS FORTRAN IV 360N-F0-'*79 3-1 


DATE 08/03/71 


TIME 01.30.50 


PAGE 0003 


C READ DATA CAROS . . _ . 

C 

00 RF.AOt l, 10021 IX, IY, IZO, IZ1 

.100? FORMAT ( 4110 ) 

IF ( IZO .SO. 0 1 GO TO 35 
I = I + l 

... : X = IX / 100. . 

Y = IY / 100. 

OX = X - X M I M 

0 Y = Y - Y M I N .. 

XPin = 1.5 * ( SORT ( DX * DX + DY * DY ) / CIN 1 
7 0 = I ? 0 - 5 3000. 

y*m n = izo / moo. t * 2.5 + 2.0 

m ( YP ( s ) .GT. 10.0 ) YPII) = 10.0 
Zl = 1 71 

Ifl 1 1 = ( Z 1 . - A M 1 N .,)_/„ D A * 3 . 5 6.0 

ip f zpm .gt. 10.0 ) zptn = 10.0 

GO TO ?0 

35 X - -1.5 +■ IP * 3. 

Y = O ♦ 0 

CALL NUMBER ( X, Y, .?1, ALT , 0.0, “1 ) 

CALL SYMBOL (-0.0, -0.0, -0.0, * FT * ,_0 . 0 , .5 » 

CALL SYMBOL (-0.0, -0.0," . 2 1 , I SYM ( I P ) , 0.0, -1 ) 

Y - 0.6 

X = -1.0 *■ IP * 2. . 

CALL SYMBOL ( X, Y, .21, NAME , 0.0, S 1 


r PLOT MAGNETIC PROFILES __ _ 

c ‘ "" ' ” . 

CALL PLOT(XPtl), Y P { 1 ) , .3) 

on 40 j = ii i _ ' 

CALL SYMBOL ( XP(J), YPtJi, .07, ISYK(IP), 0.0, -2) 

40 cn.'si r iniif 

IF ( MALT .EG. 0 ) GQ TO 2 0 

C 

C PLOT PROFILE OF ALTITUDE 

CALL PLOT ( XP( 1 > ♦ ZPtll , 3) 

00 50 J - 1, I 

CALI. S Y MB 01. ( X R f J } , Z P ( J ) , . 0 7 I S YM t I P > , .0 . 0 , . -2 . .) _ 
50 CONTINUE 
r,n T.n 20 

C • 

C AT END .O p FILE 
C 

9 00 WRITE (3,3000) ... __ 

3000 FORMAT < • END Of JOB * ) 

CALL PLOT ( 0.0, 0.0, 99 ) 

CALL EXIT . 

END 




i o - .* /v i - L 


K A f MA 1 2 


I ) A l c 


0 l / 26/ I l 


} l .-ifc 


JS 


• - T 


3 V 3 6 i 


C 

C PROGRAM TO Or T MAGNETIC DATA READY FOR SIAMPtDE PROGRAM. X AN 0 Y 
C VALUES GIVEN ONLY AT FIRST AND LAST POINT - VALUES IN BETWEEN NEED TO 
C B h CALCULATED. AVERAGES r UR EACH ALTITUDE AND OVERALL AVERAGES AN D 

C ANOMALIES ARE calculated 

C 


GOO! 


DIMENSION ' FLO (50)," 

I MF ( 5 0 ) 

0002 


01 MENS ION XI 100), 

Y(IOO) , 

0003 


REAL NAME, 

NAMEC , 


C 



0 00 A 

c 

DATA CLANK , NAMEX /' 

1 

0 0 03 


EOF = 0. 


0 0 0 6 


I =0. 


0 30 7 


00 5 J = 1,50 


0 00<i 


F L D ( J ) =0.0 


GOOD 


I M F ( J ) = 0 


U 0 1 0 


6 CUUTINUE- 



UO'l LOO) f IZ1 (100) , 1121 100) 
BLANK, OAT t , NAMEX 


EXTRA * / 


C "" 

C READ LIMITS OF DATA TO ACCEPT IF THERE ARE ANY - JUST FIELD DATA 
C A PLANK CARD WILL MEAN ALL FIELD VALUES WILL BE INCLUDED 
C 


00 I 1 
o i :i 
0013 

0 0 1 4 
GC 1 S 

00 1 o 


READ ( It 11 ) MINE, MAXF ' 

U r- i ! R M A T ( 2 1*3 ) 

UK = MA X F - HINF '• ■ 

IF ( MAXF .tO. 0 ) MA XF = 99V 99 0 
If 1 ! JX .NO. 0 ) WRITE ( 3,33) Ml Nr , MAXF 

33 FORMAT ( ' FIELD VALUES LESS THAN *, 16, * AND GREATER THAN *, 
1 16, « WILL NOT BE INCLUDED IN THE AVERAGES *) 


C 

C READ FIRST RECORD - SHOULD HAVE INITIAL X AND Y 

C 


0 0 3 7 


k F A OIL , 10 , END 

= 9 0 ) 

NAME 

? D I K t I f ’ c t XI? 

Ylt IZ» IDZ, IGE, 

DATE, 

0 0 1 :-j 

10 

FORMAT ( A 6 , 

F5 ♦ 1 

, 19, 

F6.2/ F6‘.2t IS, 

, “13, 14, • 20X, At), 

A 4 • ) ” 

0 019 

15 

NAMEC - NAME 






0»?r> 


JFE = IFF 






0021 


HD 7 = I 0 7.' 






00 >2 

60 

1 = 1 + 1 






00 3 i 


IH IFF . N t . 

0 ) 

JFE = 

IFE 



0 0 2 A 


X { I ) = XI 






0 0 25 


Y < I) = Y I 



* 



0026 


IF ( lUZ . c Q . 

0 ) 

IDZ 

= HD Z 




0027 ' 3/0(1 > = U ♦ IDZ 

0026 1/1(1) * IGE 


C . 

C TEMPORARILY SI ORE FL I GHT 'ALTITUDE IN 1 12 SPACE- NEEDED LATER WHEN 
C. A CHANGE IN FLIGHT ALTITUDE GIVES RISE TO A DIFFERENT FLIGHT NAME 
C 


0029 

' J A = JFE / 

100 ' 




00 30 
00 31 

IZ2< M = JA 
I F ( I Z 0 t I ) 

.GT. MAXF ! 

i GO 

TD 

50 

00 32 

I F ( I ZO ( I) 

.LT. M INF I 

1 GO 

TO 

50 


C ' 

C SUM FIELD FOR AVERAGING 


TIM 


00 33 
0 0 34 


F L 0 ( J A ) 
I M F ( J A ] 


FLIXJA) + IZOtl) 
1MFUAJ + I 


C READ LATER RECORDS 
C 


OATt t 


JO 3 6 

IF 

{ NAME 

.pg. blank 

) GO 

TO 

40 









00 37 

IF 

{ NAME 

.GO. NAMEC 

> GO 

TO 

40 









0 0 33 

60 S P 

= I - 

L 












DC 39 

OX 

= XU) 

- X( 1 ) 












0040 

UY 

= YII) 

- Y( 1 ) 













C STORE OAT A ON TAPE - FIRST ITEM IS THE NUMBER OF POINTS FOR 
C 


FL IGHT 


004 1 

WRITE (3) I 

CO 4? 

DO 65 J = 1,1 

00 43 

FAC = (J - 1) / 

00-4 

X ( J ) = XU ) *- F 

0 0 4 s 

Y ( J 1 = Y( I) + r 

0046 

.< L = XU) * 100 

Q0 4 7 

K2 = YU J * 1 JO 


C STORE OAT A - THE ACTUAL FLIGHT DATA 

c " 

W RITE ( A ) Kit K2 , iZO(J), I2HJ), IZ2IJ). NAMEC, J 
65 CONTINUE 

IF ( EOF .EQ. 1 . 1 GO TO 99 "" " 

I = 0 
GO TO IS 

" "90 EOF = I . """ 

GO TO 60 
99 ENDFILE 8 

REWIND 3 " ■ * 

C 

C COMPUTE AVERAGES £» LIST THEM 

X " " ' ' 

GFt.0 = 0.0 
IN DR = 0 

WRITE { 3 , 3001 ) " 

3001 FORMA T ( M • , /// 10X, ’FLIGHT ALT AVG FIELD NOR POINTS', // ) 
DO 105 J = 1 * SO 

" IF ( IMF ( J) .ED. O') GO TO 105 ”7 ” 

GELD = 'Sr t.. D +■ flbU) 

1N4R - I MAR + IMr(J) 

F LO { J ) = FLO(J) / IMF C J > 

ALT a J * 100. 

rfRITE 13,3002) ALT, rLO(J), I MF t J ) 

3002 FORMAT (• 0* , 1 2X , F6.0, 6X , F7. 0 , 7X » ~ 14 “) 

I OS CONTINUE 

GFLD - GFlD / IN OR 

WRITE! 3. 3022 ) GFLD, INBR 

3022 FORMAT { 'O' //, 12X, 1 ALL 6X, F7.0, 6X, 15 ) 

C 



ODS r Jk F S AN IV 360,\-fo-4?9 3-1 


K A T M A I 2 


REREAD DATA 
FIELD I zu 
FIELD i l 1 
FIELD ll 2 
FitLO f z 3 


07/28/71 


' 107 ( 5 , r NO = 900 } ' (V AP 

ICNT = o 

00 ISO I = l fNBR 

'<EAD 13) KL, K2 * I/O 

KA 1 j A * Nfc * 1 J G ° TG ^0 
I XTRA =0 
NAM EC = r-y A M £ 

^-i 7E <3,3003} MAMEC 
-003 FukMAT IMS' FLIGHT 
1 20 


-^1 Vfr^ CULATf: AM3MALIES' “ ~ - 

~ ACTUAL FI El. D VALUES 

I A C r T nf, L G T JUt *° t'LEVATICNS 

■ ACTUAL FlltS VALUE - Ovti'itu F0R THAT-pUGHT 
ALUt OVtkALL AVERAGE 


r - n 1 1 ! * 1 1 1 ( 1 1 * 'OA NAM E” 


110 IF £ 

>< A = 


'♦ /// ) 

KA JA } GO TO 125 


* « A6 t / / / , 1 

21 


IXTRA = fXT.RA + 

i cm r = o 

NA ' M EC = ft A M E X 

•'Fite i 3 1 3 c 0 3 ] 
125 I CUT = ICNT + 1 
is. 2 tn « i/ouj 


NAMEC 
- F'LJ(JA) 


T V _ r . y , - . 1 

.t 7 ; 4 7 I i J U ) - c F L 0 
R I 7 !.. ( j t 3 0 u 4- ) K l . K '* 1 7/1 1 

^•fmah Mio, nxl ,Zla,ri ^mriZBV NAME;- 

on ro 1 50 

FORMAT! ’ if in~ \ nv * ?? ’ . IZOfI) ’ ^ 1 I I > , 1 22 (1) 3 T NA M£ f 


J 0RMAr{ 61 l Ox] 

150 CONTINUE 

00 TO 107 

*00 Call exit 
END 


A 6 , lx, 13 ) 



